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A  method  that  connects  measurements  of  radiometric  forces  on  a  heated  vane  in  the  transitional  flow 
regime  with  the  kinetic  modeling  of  the  flow,  and  derives  the  accommodation  coefficients  through 
the  successive  analysis  of  measured  and  computed  results,  is  proposed.  The  method  utilizes  the  fact 
that  radiometric  forces  exerted  on  heated  objects  immersed  in  rarefied  gases  are  governed  by  the 
interaction  of  gas  molecules  with  the  surface.  Experimental  results  on  radiometric  forces  on  a  0. 1 1 
m  diameter  circular  vane  are  obtained  on  a  nano-Newton  thrust  stand  in  a  3  m  long  vacuum  chamber 
for  pressures  ranging  from  approximately  0.01  to  1  Pa.  The  vane  was  heated  to  419  K  on  the  hot  side 
and  396  K  on  the  cold  side.  The  numerical  modeling  is  conducted  using  a  combined  ellipsoidal 
statistical  Bhatnagar-Gross-Krook/direct  simulation  Monte  Carlo  approach  that  allows  accurate  and 
time  efficient  analysis  of  radiometric  forces  on  a  vane  in  large  vacuum  chambers  filled  with  rarefied 
gas.  Accommodation  coefficients  for  the  Maxwell  model  are  estimated  for  argon,  xenon,  and  helium 
on  a  machined  aluminum  surface,  and  found  to  be  0.81,  0.86,  and  0.53,  respectively.  ©  2009 
American  Institute  of  Physics.  [DOI:  10.1063/1.3187932] 


I.  INTRODUCTION 

The  history  of  accommodation  coefficients  of  energy  and 
momentum  of  gas  molecules  colliding  with  solid  surfaces 
spans  well  over  a  century. 1  Its  beginning  dates  to  the  work  of 
Kundt  and  Warburg"  who  studied  the  effect  of  gas  density 
change  in  the  damping  of  a  vibrating  disk.  The  viscosity 
appeared  to  decrease  with  density,  which  seemed  unexplain¬ 
able  at  the  time.  The  authors  suggested  an  incomplete  inter¬ 
action,  or  accommodation,  of  gas  molecules  at  the  surface, 
where  a  low  density  gas  slips  over  a  surface.  Following  that 
work,  Maxwell  showed  that  the  slip  phenomenon  has  roots 
in  kinetic  theory,  and  he  treated  the  solid  wall  as  something 
intermediate  between  a  perfectly  reflecting  and  a  perfectly 
absorbing  surface.  He  proposed  that  “of  every  unit  of  area  a 
portion  a  absorbs  all  the  incident  molecules,  and  afterwards 
allows  them  to  evaporate  with  velocities  corresponding  to 
those  in  still  gas  at  the  temperature  of  the  solid,  while  a 
portion  1  -a  perfectly  reflects  all  the  molecules  incident  upon 
it.”3 

The  model  proposed  by  Maxwell  is  in  fact  the  first  the¬ 
oretical  model  that  describes  gas-surface  interaction,  and  it  is 
still  widely  used  today  both  in  experiment  and  numerical 
simulation.  According  to  the  Maxwell  model,  the  velocity 
distribution  function  of  reflected  molecules  may  be  written 
as  a  function  of  the  accommodation  coefficient  a  (see,  for 
example.  Ref.  4), 

/r0,x,vr)  =  (1  -  a)fj(t,x,\r  —  2(\ r  •  n)n) 

/  02\  3/2 

+  e^2\  (1) 


denotes  the  surface  normal.  Subscripts  i  and  r  refer  to  inci- 
dent  and  reflected  molecules,  respectively,  and  /3=yml  2kTr. 
The  first  term  in  Eq.  (1)  refers  to  specular  reflection  and  the 
second  term  refers  to  diffuse  reflection.  The  reflected  tem¬ 
perature,  Tn  is  the  wall  temperature,  Tw,  according  to  the 
original  Maxwell’s  idea,  but  may  generally  be  a  free  param¬ 
eter  of  the  model. 

The  tangential  momentum  transferred  to  the  surface  by 
the  incident  molecules  may  be  written  as4 

p  ir=-m\  f,ylT(y,  ■  nMv,,  (2) 

J  vn<0 

where  the  subscript  t  refers  to  the  tangential  to  the  surface 
components  of  molecular  velocity,  and  the  tangential  mo¬ 
mentum  of  reflected  molecules  is  then 


P/  t  =  -  m  f pi rr(y r  '  n)dv,.  =  ( 1  -  a) P,v.  (3) 

J  vn>0 

These  equations  show  that  the  accommodation  coefficient 
used  in  Eq.  (1)  may  be  considered  as  the  coefficient  of  ac¬ 
commodation  of  the  tangential  momentum,  and  may  be  writ¬ 
ten  as 


(4) 


where  /  is  the  distribution  function,  t  is  time,  x  and  v  are  The  accommodation  coefficients  for  the  normal  momentum 

molecular  position  and  velocity  vectors,  respectively,  and  n  and  energy  may  be  introduced  similar  to  Eq.  (4)  as 
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a 


n 


(5) 


where  the  subscript  w  refers  to  the  surface  properties  or  the 
properties  that  would  have  had  a  gas  at  equilibrium  with  the 
wall. 

Although  the  Maxwell  model  is  still  the  most  widely 
used  model  of  gas-surface  interaction,  other  models  have 
also  been  proposed.  Here  we  mention  only  a  two-parameter 
Cercignani-Lampis  model5  that  uses  two  accommodation 
coefficients,  aT  and  and  a  multiparametric  Nocilla  model6 
in  which  the  velocity  of  reflected  molecules  is  simulated  by 
the  function 


f=nrTT  3/2cr3  exp{-  cr2[|-  crS,.]2}, 


cc  =  ' \2k!mTr , 


sr=hcr. 

The  four  parameters  Sr=  (Snr,STr),  Tn  nr  of  this  function  are 
determined  from  experimental  data. 

The  development  and  utilization  of  different  gas-surface 
interaction  models  is  related  to  various  application  areas 
where  such  interactions  are  important.  One  area  of  interest  is 
the  high  altitude  aerodynamics,  and,  in  particular,  free  mo¬ 
lecular  aerodynamics  of  satellites  (see,  for  example.  Ref.  7). 
For  the  latter  application,  the  Nocilla  model  is  often  used. 
The  importance  of  the  gas-surface  interaction  model  in  this 
case  is  obvious  since  the  collisions  of  molecules  with  the 
spacecraft  surface  are  the  dominant  process  that  influences 
drag,  lift,  and  heat  loads. 

Another  area  where  the  gas-surface  processes  are  impor¬ 
tant  is  gas  flows  in  micro-  and  nanoscale  devices.  In  such 
devices,  the  gas  mean  free  path  is  comparable  to  character¬ 
istic  flow  dimensions,  and  the  consideration  of  kinetic  effects 
is  essential  for  accurate  prediction  of  device  performance  and 
peculiarities.  The  large  surface-to-volume  ratio  further  in¬ 
creases  the  influence  of  the  wall.  Note  that  for  microscale 
flows,  the  preservation  of  the  detailed  balance  in  collisions  of 
gas  molecules  with  solid  interfaces  is  critical.  Therefore,  the 
Nocilla  model,  which  does  not  satisfy  this  requirement,  is  not 
a  good  choice,  and  the  Maxwell  and  Cercignani-Lampis 
models  are  better  suited  for  the  description  of  low  speed 
flows  in  microdevices.  Beside  these  two  areas,  gas-surface 
interaction  is  important,  if  not  determining,  in  many  other 
applications.  Near-continuum  supersonic  flows  over  sharp 
leading  edges,  contamination  problems,  and  two-phase 

o 

flows1  are  just  a  few  examples  of  such  applications. 

Accurate  prediction  of  the  above  flows  requires  the  re¬ 
searcher  not  only  to  select  an  appropriate  gas-surface  inter¬ 
action  model,  but  also  to  specify  the  parameters  of  this 
model  for  each  type  of  gas  species-solid  wall  interface.  Two 
principal  approaches  are  used  to  determine  parameters  of  the 
model,  theoretical  and  experimental.  The  theoretical  ap¬ 
proach  is  usually  based  on  the  detailed  studies  of  molecular 
interactions  using  classical  or  quasiclassical  trajectory  calcu¬ 
lations  in  the  framework  of  the  molecular  dynamics  method.9 
In  the  experimental  approach,  parameters  of  the  selected  in¬ 


teraction  model  are  estimated  directly  from  the  measure¬ 
ment.  The  parameters  for  the  Nocilla  model,  for  example,  are 
usually  obtained  from  molecular  beam  experiments  or  flight 
experiments  (see,  for  example.  Refs.  10  and  11  and  refer¬ 
ences  therein).  In  Ref.  12,  a  connection  between  the  exit 
velocity  distribution  described  by  the  Nocilla  model  and  the 
classical  momentum  and  energy  accommodation  coefficients 
was  given. 

The  advantage  of  the  molecular  beam  technique  is  that  it 
may  provide  detailed  information  on  the  velocity  distribu¬ 
tions  of  reflected  molecules.  There  are  many  situations,  how¬ 
ever,  when  such  detailed  information  is  not  necessary,  and 
the  knowledge  of  accommodation  coefficients,  either  mo¬ 
mentum  or  energy,  would  suffice.  Examples  include  the  force 
estimate  of  spacecraft  at  high  altitudes,  or  the  evaluation  of 
heat  loads  in  microdevices.  Over  the  past  three  decades,  mo¬ 
lecular  beam  experiments  have  been  used  extensively  to  de¬ 
termine  both  the  normal  and  tangential  momentum  and  en¬ 
ergy  accommodation  coefficients 13-1 5  for  various  gas-surface 
pairs.  For  the  energy  (thermal)  accommodation  coefficient, 
parallel  plates,  coaxial  cylinders,  and  hot-wire  methods  have 
been  widely  used.  A  comprehensive  review  of  different  ap¬ 
proaches  to  the  thermal  accommodation  coefficient  measure¬ 
ments  may  be  found  in  Ref.  16.  A  wide  range  of  results  for 
accommodation  coefficients  have  been  reported  for  the  three 
gases  studied  here  (argon,  xenon,  and  helium).  In  view  of 
this,  the  results  of  Ref.  17  are  particularly  important;  these 
are  state-of-the-art  measurements  of  energy  accommodation 
coefficients  for  argon  and  helium  at  temperatures  similar  to 
those  for  the  present  work,  and  thus  provide  the  most  useful 
comparison  with  the  present  results.  Various  experimental 
techniques  used  in  the  past  to  measure  tangential  momentum 
accommodation  coefficient,  such  as  the  rotating  cylinder 
method,  the  spinning  rotor  gage  method,  the  flow  through 
microchannel  approach,  as  well  as  the  molecular  beam  tech¬ 
nique,  are  discussed  in  recent  review  article.18 

In  contrast  to  high-enthalpy  flows  around  space  vehicles, 
gas-driven  flows  in  microscale  devices  are  characterized  by 
relatively  low  gradients  in  gas  velocity  and  temperature,  and 
the  velocity  distribution  function  in  these  flows  is  often  close 
to  Maxwellian.  As  a  result,  prediction  of  gas-driven  flows  in 
such  devices  typically  requires  knowledge  of  momentum 
and/or  energy  accommodation  coefficients  as  a  function  of 
gas  and  surface  temperature.  The  use  of  molecular  beam 
technique  may  be  quite  difficult  in  this  case  since  the  after¬ 
collision  velocities  need  to  be  analyzed  for  a  large  number  of 
precollisional  energies.  On  the  other  hand,  standard  tech¬ 
niques  for  the  accommodation  measurement  may  not  be  ap¬ 
plicable  when  information  on  momentum  accommodation, 
normal,  or  tangential  is  needed. 

One  major  issue  with  the  need  of  gas-surface  interaction 
parameters  to  predict  complicated  flow  interactions  is  the 
range  of  experimental  data  for  similar  flows.  For  example, 
various  experiments  can  have  different,  and  in  some  cases 
conflicting,  results.  Take,  for  instance,  the  measurement  of 
the  energy  accommodation  of  helium  on  a  platinum  surface. 
References  19-21  present  results  which  vary  by  more  than 
30%.  Thus  there  is  a  need  to  reinvestigate  these  data  sets  in 
order  to  study  the  effects  of  gas  temperature,  surface  tem- 
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perature,  surface  preparation,  gas  adsorption  on  surfaces,  and 
gas  pressure. 

The  objective  of  this  work  is  to  evaluate  the  feasibility  of 
using  a  new  method  for  measurements  of  momentum  accom¬ 
modation  coefficients,  based  on  the  comparing  experimental 
and  computational  results  on  radiometric  forces  on  heated 
plates.  Radiometric  forces  are  typically  exerted  on  nonuni- 
forrnly  heated  objects  immersed  in  rarefied  gases,  and  tend  to 
move  these  objects  in  the  direction  from  the  cold  to  the  hot 
side.  The  authors  of  Ref.  22  have  recently  remarked  that  the 
measurement  of  radiometric  forces  may  yield  data  on  gas- 
surface  interaction.  However,  there  are  two  problems  that 
make  the  direct  use  of  such  measurements  to  infer  momen¬ 
tum  accommodation  coefficients  extremely  difficult.  First, 
there  are  usually  molecular  collisions  present  in  radiometric 
flows,  and  these  collisions  do  not  allow  simple  and  accurate 
analytic  evaluation  of  accommodation  coefficients  from 
force  measurements  beyond  the  free  molecular  regime.  Sec¬ 
ond,  while  the  availability  of  force  measurements  in  free 
molecular  regime  would  offer  the  benefit  of  accommodation 
coefficient  evaluation,  there  is  a  physical  limitation  in  the 
accuracy  of  such  measurements.  The  fewer  gas-surface  col¬ 
lisions  a  radiometer  vane  experiences,  the  greater  the  experi¬ 
mental  error.  To  avoid  these  difficulties,  it  is  suggested  in  the 
present  work  to  measure  radiometric  forces  in  the  transitional 
regime,  and  then  use  kinetic  modeling  of  radiometric  flows 
to  infer  the  momentum  accommodation  coefficients. 


II.  RADIOMETRIC  APPROACH  TO  MOMENTUM 
ACCOMMODATION  STUDY 


fhJyr  ■  n)rfv,.  =  n 


I 


fjyr  ■  n )d\r 


(7) 


where  subscript  g  refers  to  the  incident  gas  molecules.  In  a 
free  molecular  flow,  the  Maxwellian  distribution  function  of 
fg  may  be  reasonably  assumed.  For  fully  diffuse  accommo¬ 
dation,  the  number  density  of  reflected  molecules  is  obtained 
by  integrating  Eq.  (7)  over  equilibrium  distribution  functions 
to  give 


where  indices  h  and  c  refer  to  either  hot  or  cold  side  of  the 
plate.  The  force  on  the  side  of  the  plate  will  be 


where  pg  is  the  gas  pressure.  The  first  term  in  Eq.  (9)  is  for 
the  incident  molecules  and  the  second  term  accounts  for  the 
contribution  from  the  reflected  molecules.  If  the  accommo¬ 
dation  coefficient  is  introduced  according  to  the  expression 

23 

suggested  by  Knudsen,” 


(10) 


then  assuming  the  same  accommodation  coefficient  on  the 
hot  and  cold  sides  of  the  plate  (a  small  temperature  differ¬ 
ence  between  the  plates),  and  using  Tw  from  Eq.  (10)  instead 
of  Th  and  T,  in  Eq.  (9),  one  can  obtain  the  expression  for  the 
total  radiometric  force  on  the  plate. 


The  radiometric  forces  on  a  heated  plate  may  be  de¬ 
scribed  analytically  only  in  a  free  molecular  regime;  the 
presence  of  even  relatively  small  number  of  molecular  colli¬ 
sions  in  the  transition  regime  complicates  the  flow  to  the 
point  where  accurate  analytical  description  is  not  possible, 
and  a  numerical  approach  has  to  be  used  to  address  the  prob¬ 
lem.  Even  for  a  free  molecule  flow,  some  model  needs  to  be 
used  for  the  gas-surface  accommodation  in  order  to  make 
analytical  treatment  possible. 

Generally,  for  a  plate  with  its  opposite  sides  heated 
uniformly  to  different  temperatures  Th  and  Tc,  the  forces  in 
the  direction  normal  to  the  plate,  created  by  molecules  re¬ 
flected  from  the  hot  and  the  cold  sides  of  the  plate,  may  be 
written  as 


r  pJ  l(l~  aK)Tg+  aKTh 
Fk~  2\\  Th 

_  ^/(  1  ~  aK)Tg  +  QifTc | 


(ID 


Note  that  the  contributions  from  the  incident  molecules  can¬ 
cel  out  in  the  free  molecular  flow;  the  force  is  directed  from 
the  hot  to  the  cold  surface. 

If  the  Maxwell  model  of  gas-surface  interaction  is  used, 
then,  substituting  Eq.  (1)  into  Eq.  (6)  and  making  use  of  Eq. 
(7),  one  can  obtain  for  the  free  molecular  force 


Fh  =  nhm 


Fc  =  ncm 


ftyjyr '  n  )d\r 

/cvrT(vr  •  n  )d\n 


and 


(6) 


where  subscripts  h  and  c  refer  to  the  hot  and  cold  sides, 
respectively.  The  number  density  that  describes  the  flux  of 
reflected  molecules  may  be  obtained  from  the  assumption  of 
the  equality  of  the  incident  and  reflected  mass  flux  (i.e.,  no 
sticking  on  the  surface). 


Therefore,  the  free  molecular  radiometric  force  calculated 
using  the  Maxwell  model  is  linearly  dependent  on  the  tan¬ 
gential  momentum  accommodation  coefficient.  If  the  mo¬ 
mentum  accommodation  coefficient  in  the  Maxwell  model, 
a,  is  close  to  the  energy  accommodation  coefficient  in  the 
Knudsen  model,  aK,  then  the  force  predictions  obtained  with 
Eqs.  (11)  and  (12)  are  similar  for  small  temperature  differ¬ 
ences.  The  difference  between  them  becomes  significant 
when  the  surface  temperatures  are  not  similar.  Equation  (12) 
allows  one  to  easily  calculate  the  accommodation  coeffi¬ 
cients  when  the  radiometric  force  in  the  free  molecular  re¬ 
gime  can  be  measured. 


Downloaded  02  Mar  2010  to  198.14.54.253.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://pof.aip.org/pof/copyright.jsp 


073101-4  Selden  et  al. 


Phys.  Fluids  21,  073101  (2009) 


In  reality,  however,  it  is  difficult  to  accurately  measure 
the  radiometric  force  in  free  molecular  regime.  Such  mea¬ 
surements  are  possible  for  a  transitional  regime,  for  which 
the  above  analytic  expressions  are  not  applicable.  Therefore, 
it  is  reasonable  to  infer  the  accommodation  coefficients  from 
a  numerical  simulation  performed  for  a  given  gas-surface 
interaction  model  with  varying  parameters  of  the  model.  It  is 
clear  that  the  conventional  continuum  approaches  of  the 
computational  fluid  dynamics,  such  as  those  based  on  the 
solution  of  the  full  Navier-Stokes  equations  or  boundary 
layer  equations,  cannot  be  used  to  compute  radiometric 
forces  in  the  transitional  flow  regime.  In  these  approaches, 
developed  for  modeling  gas  flows  close  to  equilibrium,  the 
effects  of  rarefaction  are  typically  accounted  for  through  the 
boundary  conditions  of  slip  velocity  and  temperature  jump 
on  the  surface.  The  assumption  of  small  deviation  from  equi¬ 
librium  makes  them  inapplicable  for  modeling  radiometric 
flows  and  calculating  radiometric  forces.  In  this  case,  a  ki¬ 
netic  approach  based  on  the  solution  of  the  Boltzmann  equa¬ 
tion  has  to  be  used.  For  a  kinetic  approach,  a  kinetic  model 
of  gas-surface  interaction  needs  to  be  used,  such  as  the  Max¬ 
well  model,  and  the  approach  naturally  gives  the  velocity 
distribution  functions  for  the  incident  and  reflected 
molecules. 

Thus,  in  order  to  obtain  the  accommodation  coefficients 
for  a  given  gas-surface  interaction  model,  numerical  results 
need  to  be  obtained  by  comparing  results  from  a  kinetic  ap¬ 
proach  for  a  computational  setup  that  closely  reproduces  the 
experimental  one.  The  details  on  the  present  experimental 
setup  are  given  in  Sec.  III. 

III.  EXPERIMENTAL  SETUP 

As  radiometric  phenomena  occur  in  rarefied  conditions, 
there  are  only  two  ways  to  study  them  experimentally.  The 
first  is  to  build  extremely  small  devices  on  the  order  of  na¬ 
nometers  and  test  them  under  atmospheric  conditions.  The 
second  way  is  to  build  a  larger  device  and  modify  the  back¬ 
ground  pressure  such  that  the  local  Knudsen  number  is  large 
enough  for  the  flow  to  be  considered  transitional  (i.e.,  Kn 
>0.01).  In  this  work  the  latter  method  has  been  chosen,  and 
all  the  experimental  results  that  follow  have  been  achieved 
under  low  pressure  conditions  in  a  large  3.0  m  diameter 
vacuum  chamber.  The  use  of  such  a  large  chamber  is  critical 
to  avoid  the  effect  of  chamber  walls  that  was  found  to 
strongly  impact  the  radiometric  force  in  smaller  chambers."4 

To  accurately  measure  the  impact  of  various  accommo¬ 
dation  coefficients,  and  to  be  practical  to  model  using  an 
axisymmetric  code,  a  circular  radiometer  vane  with  a  diam¬ 
eter  of  11.13  cm  was  used.  The  vane  consisted  of  a  Teflon 
insulator  sandwiched  between  two  aluminum  plates  with  a 
resistive  heater  located  between  one  of  the  plates  and  the 
insulator.  The  temperature  of  one  side  of  the  device  was 
maintained  by  varying  the  power  input  to  the  heater,  while 
the  temperature  of  the  opposite  side  was  not  actively  main¬ 
tained  and  was  allowed  to  float.  Each  of  the  three  pieces  of 
the  radiometer  vane  had  a  thickness  of  0.32  cm,  and  when 
assembled  yield  a  total  device  thickness  of  0.96  cm. 

One  motivation  for  this  particular  configuration  of  radi- 


FIG.  1.  Setup  of  the  radiometric  experiment. 


25 

ometer  vane  comes  from  historical  work"  where  rudimen¬ 
tary  temperature  measurements  of  the  vanes  suggested  that  a 
significant  temperature  drop  occurred  at  the  outermost  edges. 
This  same  work  made  it  quite  clear  that  to  accurately  deduce 
a  theory  for  the  operation  of  the  radiometer,  it  would  be 
necessary  to  discover  exactly  what  effect  the  temperature 
variations  at  the  edges  had.  For  the  sake  of  clarity  it  should 
be  noted  here  that  there  are  two  gradients  important  to  the 
flow:  The  first  of  these  shall  be  referred  to  as  the  radial 
gradient  and  will  refer  to  the  temperature  profile  of  a  plate 
from  the  center  to  the  periphery,  while  the  second  will  be 
called  the  axial  gradient  and  will  refer  to  temperature  profile 
along  an  axis  normal  to  the  face.  In  an  ideal  experiment,  the 
axial  gradient  would  be  large  and  the  radial  gradient  would 
be  nonexistent  such  that  the  experiment  and  simulation  share 
nearly  identical  temperature  profiles.  It  is  for  these  reasons 
that  the  particular  aluminum  “sandwich”  design  was  chosen; 
not  only  does  the  high  thermal  conductivity  maximize  the 
surface  temperature  of  the  hot  plate  (and  thus  the  axial  tem¬ 
perature  gradient),  but  it  also  minimizes  the  radial  tempera¬ 
ture  gradients  near  the  edges  of  the  device. 

To  measure  the  force  produced  by  this  device,  it  was 
mounted  on  a  modified  nano-Newton  Thrust  Stand"'1  (nNTS) 
located  inside  the  vacuum  chamber.  A  schematics  of  the  ex¬ 
perimental  setup  is  shown  in  Fig.  1.  Here,  every  effort  was 
made  to  minimize  the  impact  of  the  thrust  stand  arm  and 
attachment  mechanism  by  using  6.35  mm  tubing  coupled 
with  a  2  X  40  mm2  threaded  rod.  When  calibrated  using  a 
pair  of  electrostatic  combs,”  the  nNTS  provides  very  accu¬ 
rate  and  repeatable  data  with  typical  force  resolution  of  ap¬ 
proximately  0.1  /aN  and  statistical  scatter  of  about  1%.  For 
the  preliminary  experiment,  the  experimental  error  based  on 
standard  deviation  ranges  from  a  few  percent  at  the  lowest 
pressures  to  less  than  1%  through  most  of  the  curve.  How¬ 
ever,  due  to  the  normalization  by  experimental  temperature 
measurements  and  the  small  uncertainty  of  the  calibration 
method,  the  total  absolute  experimental  uncertainty  is  —4%. 
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Day-to-day  variation  of  multiple  data  sets  has  been  observed 
to  be  —1%. 

The  experimental  data  were  obtained  by  evacuating  the 
vacuum  chamber  to  a  base  pressure  below  10-3  Pa.  This  low 
pressure  was  required  to  minimize  the  impact  of  the  back¬ 
ground  gas  to  a  level  low  enough  as  to  be  inconsequential  to 
the  measurements  being  made.  While  the  evacuation  of  the 
chamber  was  taking  place,  a  constant  voltage  was  applied  to 
the  heater.  This  resulted  in  the  main  radiometer  surfaces 
reaching  temperatures  of  approximately  419  K  (hot)  and  394 
K  (cold),  although  the  exact  values  fluctuated  depending  on 
both  the  species  and  pressure  of  the  background  gas.  Force 
measurements  were  made  by  varying  the  background  pres¬ 
sure  of  the  gas  in  the  chamber,  where  argon,  helium,  and 
xenon  were  all  used.  The  highest  background  pressure 
achieved  was  approximately  1 .6  Pa,  but  varied  depending  on 
the  molecular  weight  of  the  background  gas. 


IV.  NUMERICAL  MODELING  OF  RADIOMETRIC 
FLOWS:  A  COMBINED  KINETIC  APPROACH 

In  this  work,  a  combined  ellipsoidal  statistical 
Bhatnagar-Gross-Krook/direct  simulation  Monte  Carlo 
(ES-BGK/DSMC)  approach,  where  the  final  solution  is  ob¬ 
tained  in  two  successive  steps.  First,  an  ES-BGK  modeling  is 
conducted  in  a  large  computational  domain  that  includes 
both  the  radiometer  vane  and  the  chamber  walls.  The  solu¬ 
tion  of  this  first  step  is  used  to  set  the  boundary  conditions 
for  the  second  step.  At  the  second  step,  the  DSMC  method  is 
applied  in  a  much  smaller  domain,  with  the  subsonic  bound¬ 
ary  conditions  taken  from  the  first  step.  The  use  of  such  a 
new  approach  is  based  on  the  fact  that  DSMC  modeling  of  a 
radiometric  flow  on  a  10  cm  vane  in  a  3  m  chamber,  where 
the  accuracy  of  the  radiometric  force  modeling  needs  to  be 
on  the  order  of  1%,  is  prohibitively  expensive  even  for  mod¬ 
ern  parallel  computers.  On  the  other  hand,  the  ES-BGK 
method  was  found  to  be  fairly  accurate  in  predicting  all  gas 
macroparameters  in  the  computational  domain,  but  overpre¬ 
dicting  the  DSMC  results  on  radiometric  forces  by  ~  10%  in 
the  range  of  pressures  where  the  force  is  near  its  maximum. 
This  is  related  to  the  approximations  inherent  in  the  ES-BGK 
equation,  and  difficulty  of  modeling  the  radiometric  force, 
that  is  typically  less  than  1%  of  the  force  on  either  cold  or 
hot  side  of  the  vane. 

In  this  work,  the  computational  tool  SMILE  (Ref.  28)  was 
used  to  obtain  the  solutions  with  the  DSMC  method.  In 
DSMC  runs,  the  variable  soft  sphere  model  with  parameters 
listed  in  Ref.  29  was  used  for  the  molecular  collisions,  and 
the  Maxwell  model  was  used  to  calculate  gas-surface  colli¬ 
sions.  A  finite  volume  solver  SMOKE  (Ref.  30)  has  been  used 
to  deterministically  solve  the  ES  model  kinetic  equation. 
SMOKE  is  a  parallel  code  based  on  conservative  numerical 
schemes  developed  by  Mieussens.  A  second  order  spatial 
discretization  was  used.  The  solutions  were  typically  ob¬ 
tained  in  two  successive  steps.  First,  an  implicit  time  integra¬ 
tion  scheme  was  run  until  the  result  is  converged.  Second,  a 
conservative  explicit  time  integration  scheme  was  used  with 
the  initial  conditions  from  the  first  step.  This  two-step 


approach  allowed  up  to  two  orders  of  magnitude  reduction  in 
computational  time  compared  to  an  explicit-only  case. 

The  four  macroparameters  (density,  temperature,  and 
two  velocities)  from  the  ES-BGK  solution  were  used  at  the 
external  boundaries  of  the  DSMC  computational  domain. 
That  means  that  the  velocities  of  molecules  entering  the 
DSMC  computational  domain  are  sampled  from  the  Max¬ 
wellian  distribution  with  parameters  from  the  ES-BGK  solu¬ 
tion.  It  is  important  that  the  ES-BGK  macroparameters  used 
in  the  DSMC  boundary  conditions  were  computed  from  the 
incoming  fluxes  only.  Good  agreement  between  the  full 
DSMC  and  the  combined  kinetic  approach,  obtained  by  the 
authors  on  a  smaller,  0.2  m  domain,  allowed  the  application 
of  the  combined  approach  to  analyze  radiometric  flows  in  a 
large  vacuum  chamber.  A  3  m  cylindrical  chamber  is  simu¬ 
lated  in  this  work,  whose  geometry  with  good  accuracy  re¬ 
produces  the  companion  experimental  setup.  The  radiometer 
size  and  location  inside  the  chamber,  as  well  as  the  tempera¬ 
ture  conditions,  also  correspond  to  those  used  in  the  experi¬ 
ment.  Diffuse  reflection  with  a  complete  energy  and  momen¬ 
tum  accommodation  was  assumed  on  the  chamber  walls  and 
the  surface  of  the  vane  (with  one  exception  explained  be¬ 
low).  Since  the  experimental  setup  closely  approximates  a 
flow  with  an  axial  symmetry,  axisymmetric  ES-BGK  and 
DSMC  codes  were  used  in  these  computations.  The  subsonic 
boundaries  of  the  DSMC  computational  domain  were  located 
30  cm  from  the  vane  both  in  the  axial  and  radial  directions. 

V.  EVALUATION  OF  ACCOMMODATION 
COEFFICIENTS 

Three  gases  were  considered  in  this  work,  argon,  xenon, 
and  helium.  The  radiometric  forces  for  these  gases,  obtained 
with  the  combined  ES-BGK/DSMC  approach  as  well  as 
measured  experimentally,  are  presented  in  Fig.  2.  Generally, 
the  radiometric  force  consists  of  two  components,  (i)  the 
total  radiometric  force  that  includes  the  force  resulting  from 
the  pressure  difference  between  the  hot  and  the  cold  sides  of 
the  vane,  and  (ii)  the  shear  force  on  the  lateral  (circumferen¬ 
tial)  side  of  the  vane.  To  show  separate  contribution  of  these 
forces,  two  sets  of  numerical  results  are  shown,  the  total 
radiometric  force  that  includes  both  component,  and  the  ra¬ 
diometric  force  that  is  based  on  pressure  alone.  The  results 
show  that  the  shear  force  is  a  minor  factor  for  pressures 
smaller  than  0.6  Pa  for  argon,  where  the  maximum  force  is 
observed.  It  becomes  more  significant  for  larger  pressures, 
for  which  the  contribution  of  the  lateral  side  of  the  vane 
cannot  be  ignored.  It  may  appear  preferable  to  analyze  the 
accommodation  coefficients  under  conditions  where  the  lat¬ 
eral  side  has  a  negligible  effect,  such  as  a  much  thinner  vane, 
but  it  is  difficult  to  realize  in  the  experiment. 

For  all  three  gases,  the  experimental  data  lay  lower  than 
the  numerical  points,  which  is  a  clear  indication  of  an  incom¬ 
plete  surface  accommodation.  Beyond  that,  several  other  fac¬ 
tors  may  play  a  role  in  this  difference.  First,  there  are  nu¬ 
merical  and  experimental  errors;  they  are  not  expected  to 
cause  a  difference  between  the  computation  and  the  measure¬ 
ment  larger  than  5%.  Then,  there  is  a  finite  chamber  size, 
with  unknown  accommodation  on  chamber  walls.  This  has 
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FIG.  2.  Experimental  and  computed  radiometric  force  for  argon  (left),  xenon  (center),  and  helium  (right). 


been  found  to  be  a  minor  issue  in  a  series  of  ES-BGK  cal¬ 
culations,  where  the  chamber  size  larger  than  about  2  m  was 
found  to  have  a  negligible  effect  on  the  radiometric  force. 
Finally,  there  is  some  impact  of  intermolecular  collision  law, 
or,  in  other  words,  gas  viscosity  and  heat  conductivity.  This 
factor  is  also  believed  to  be  minor,  as  the  bulk  gas  properties 
correspond  to  well  established  experimental  values  for  the 
temperature  range  under  consideration.  All  these  indicate 
that  the  gas  accommodation  on  the  vane  surface  is  the  main 
reason  for  the  difference  between  the  numerical  and  experi¬ 
mental  values.  For  comparison,  additional  computations 
were  conducted  for  the  thee  gases  using  the  Maxwell  model 
with  an  accommodation  coefficients  of  0.8  for  argon,  0.86 
for  xenon,  and  0.5  for  helium.  It  is  clearly  seen  that  the  use 
of  a  lower  accommodation  coefficient  allows  one  to  obtain 
good  agreement  with  experimental  data  for  all  pressures  and 
gases  considered. 

The  value  of  0.5  for  accommodation  coefficient  used  for 
helium  to  reproduce  the  experimental  data  is  in  fact  close  to 
the  experimental-to-computed  ratio  of  0.53  obtained  after  av¬ 
eraging  over  pressures.  Remember  that  the  Maxwell  model  is 
characterized  by  a  linear  dependence  between  the  force  and 
the  accommodation  coefficient;  obviously,  the  dependence  is 
close  to  linear  in  the  transitional  regime  as  well.  The  value  of 
0.5  may  in  fact  be  obtained  if  the  Knudsen  model  of  accom¬ 
modation  is  assumed,  and  the  functional  dependence  of  Eq. 
(11)  is  applied.  In  this  case,  the  unknown  a  is  calculated  by 
equating  the  ratio  of  the  right  hand  sides  of  Eq.  (11)  with 
aK=  a  and  aK=  1  to  the  experimental-to-computed  force  ra¬ 
tio.  Note  that  the  value  0.5  is  obtained  when  the  Knudsen 
model  of  accommodation  is  assumed,  and  the  functional  de¬ 
pendence  of  Eq.  (11)  and  equating  the  ratio  of  the  right  hand 
sides  of  Eq.  (11)  with  an  unknown  aK  and  aK=  1  to  the 
experimental-to-computed  force  ratio. 

Interestingly,  the  value  of  0.5  also  coincides  with  that  of 
the  Knudsen  model  of  accommodation  obtained  assuming 
the  functional  dependence  of  Eq.  (11)  and  equating  the  above 
ratio.  As  was  mentioned  earlier,  the  difference  between  the 
accommodation  coefficients  defined  by  the  Maxwell  model 
and  the  Knudsen  expression  is  small  for  relatively  small  tem¬ 
perature  differences  examined  in  this  work.  It  is  therefore 
impossible  to  state  which  one  is  a  better  approximation  for 
the  transitional  regime.  For  kinetic  approaches,  the  authors 
believe  that  the  use  of  Eq.  (12)  may  be  a  better  fit,  with  a 


simple  ratio  between  the  experimental  and  numerical  radio- 
metric  forces  being  an  estimate  of  the  accommodation  coef¬ 
ficient  in  the  Maxwell  model.  Such  a  ratio  for  different  gases 
is  presented  in  Fig.  3.  The  accommodation  coefficients  for 
the  Maxwell  model,  obtained  in  this  work,  are  0.81  for  ar¬ 
gon,  0.86  for  xenon,  and  0.53  for  helium,  all  of  them  on  a 
machined  aluminum  surface.  Note  that  the  value  of  the  ac¬ 
commodation  coefficient  increases  with  molecular  mass, 
which  is  consistent  with  the  experimental  observation  of  Ref. 
32  but  contradicts  to  a  hypothesis  of  Ref.  33. 

VI.  COMPARISON  WITH  PREVIOUS  MEASUREMENTS 

Comparison  of  the  above  accommodation  coefficients 
with  those  measured  in  the  past  is  complicated  by  several 
factors  in  addition  to  their  obvious  dependence  on  particular 
gas  and  surface  material.  First,  the  coefficients  obtained  in 
this  work  are  integral  and  not  incident  angle  dependent. 
Therefore,  it  is  difficult  to  compare  them  with  molecular 
beam  experiments.  Second,  the  coefficients  are  generally 
sensitive  to  the  wall  and  the  surrounding  gas  temperature, 
and  the  results  should  be  analyzed  for  the  same  temperature 
regime.  Finally,  the  purity  of  the  surface  is  very  important,  as 
the  surface  coverage  and  surface  contamination  change  the 
accommodation  coefficients.  The  last  factor  is  related  to  the 
surface  temperature  and  associated  gas  desorption,  the  sur- 
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FIG.  3.  Experimental-to-computed  radiometric  force  ratio. 
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TABLE  I.  Comparison  of  the  present  accommodation  coefficients  to  the  previous  measurements. 


Gas 

Present 

Ref.  17 

Ref.  32 

Ref.  16 

Ref.  18 

Helium 

0.53 

0.38-0.47 

0.65 

Argon 

0.81 

0.86 

0.334-0.75 

0.893 

Xenon 

0.86 

0.4 

0.95 

face  roughness  and  multiple  gas-surface  encounters,  and  the 
gas  pressure.  All  these  factors  contribute  to  large  differences 
often  observed  between  accommodation  coefficients  mea¬ 
sured  by  different  authors.  Some  of  the  published  results  for 
the  three  gases  considered  in  this  work,  along  with  the 
present  data,  are  given  in  Table  I. 

For  helium,  the  present  accommodation  coefficients, 
0.53  for  the  Maxwell  model  and  0.5  for  the  Knudsen  model, 
are  close  to  that  obtained  in  Ref.  17  for  the  thermal  accom¬ 
modation  coefficient  on  machined  aluminum  kept  at  room 
temperature,  for  which  the  value  of  0.47  was  measured.  The 
accommodation  coefficient  of  helium  on  a  plasma  treated 
surface  obtained  in  Ref.  17  is  lower,  0.38.  The  normal  mo- 
mentum  coefficients  recommended  “  for  helium  on  alumi¬ 
num  are  somewhat  higher,  0.65.  It  was  also  shown  in  Ref.  32 
that  the  efficiency  of  the  momentum  transfer  process  in¬ 
creases  with  the  mass  of  gas  molecules,  and  relatively 
weakly  depends  on  the  surface  material  for  temperatures 
ranging  from  25  to  550  °C.  The  tangential  momentum  coef¬ 
ficients  of  helium  on  aluminum  are  not  available,  but  for 
other  materials  were  found  to  vary  in  a  wide  range  depend¬ 
ing  on  the  experimental  technique  used,  from  0.2  (Ref.  34)  to 
about  0.9.35 

The  thermal  accommodation  coefficient  of  argon  on  alu¬ 
minum,  tabulated  in  Ref.  16,  ranges  from  0.334  to  0.75  for 
different  experimental  techniques  and  surface  temperatures 
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from  400  to  800  K.  A  larger  value  of  0.86  was  measured 
for  argon  atoms  colliding  with  a  machined  aluminum  sur¬ 
face.  A  tangential  momentum  accommodation  coefficient  of 
0.893  was  recommended  in  Ref.  18  based  on  the  analysis  of 
a  large  array  of  experimental  data. 

The  accommodation  of  xenon  on  aluminum  has  not  been 
extensively  studied  in  the  past.  The  thermal  accommodation 
coefficient  was  reported  for  temperatures  from  500  to  800  K 
as  0.4, 16  where  a  concentric-cylinder  method  was  used. 
Among  other  materials,  platinum  was  studied  theoretically  ’6 
and  the  energy  and  momentum  accommodation  coefficients 
were  calculated  for  room  temperature  conditions  to  be  0.85 
and  0.81,  respectively.  A  mean  value  of  0.95  was  recom¬ 
mended  in  Ref.  18  for  the  tangential  momentum  accommo¬ 
dation  coefficient  of  xenon  on  commonly  employed  surface 
materials.  In  measurement/  this  coefficient  was  estimated 
as  0.9  for  xenon  on  bronze  ribbon. 

VII.  CONCLUSIONS 

A  method  for  estimation  of  gas-surface  accommodation 
coefficients,  based  on  comparing  measured  and  computed 
radiometric  forces  on  heated  vanes  in  rarefied  flows,  is  pre¬ 
sented.  The  method  applies  a  new  combined  ES-BGK  kinetic 
approach  to  match  accurately  measured  force  on  a  circular 


radiometer  installed  on  a  nNTS  and  mounted  in  a  large 
vacuum  chamber.  Accommodation  coefficients  for  the  Max¬ 
well  model  of  gas-surface  interaction  may  be  deduced  for  a 
given  pressure  and  gas-surface  pair  whether  through  the  suc¬ 
cessive  use  of  the  combined  approach  with  different  values 
of  the  accommodation  coefficient,  or  assuming  a  linear  de¬ 
pendence  of  radiometric  force  on  the  accommodation  coeffi¬ 
cient. 

Helium,  argon,  and  xenon  were  considered  in  this  work, 
for  pressures  ranging  from  approximately  0.01  to  1  Pa,  and 
an  aluminum  vane  with  a  diameter  of  0.113  m  was  exam¬ 
ined.  The  suggested  values  of  the  Maxwell  model  accommo¬ 
dation  coefficients  are  0.81  for  argon,  0.86  for  xenon,  and 
0.53  for  helium,  which  reasonably  agree  with  momentum 
and  energy  accommodation  coefficients  proposed  in  litera¬ 
ture.  The  proposed  experimental-computational  method  is 
general  enough  to  be  applied  to  a  wide  range  of  gases,  sur¬ 
faces,  and  temperature  conditions. 
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